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Abstract 

We address the construction of stable random matrix ensembles as 
the generalization of the stable random variables (Levy distributions). 
With a simple method we derive the Cauchy case, which is known 
to have remarkable properties. These properties allow for such an 
intuitive method -that relies on taking traces- to hold. Approximate 
but general results regarding the other distributions are derived as 
well. Some of the special properties of these ensembles are evidenced 
by showing partial failure of mean-field approaches. To conclude, we 
compute the confining potential that gives a Gaussian density of states 
in the limit of large matrices. The result is an hyper geometric function, 
in contrast with the simplicity of the Cauchy case. 

1 Introduction 

The ensembles of random matrices have shown beyond doubt its usefulness 
in many physical applications 13J. The definition of the joint probability 
distribution P (M) of the matrix elements of a iV by N matrix M is |13j : 

P(M) =C N exp[-Tr(V(M))] , (1) 

with an arbitrary V (M) , provided existence of the partition function Cjf. 
Integrating over the parameters related to the eigenvectors, one obtains 
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the well-known joint probability distribution of the eigenvalues: 

N 

P(x u ...,x N ) = C N exp[-^2V(xi)]Y[\xi - Xjf . (2) 

i=l i<j 

Level repulsion described by the Vandermonde determinant is originated 
from the Jacobian, that appears when passing from the integration over 
independent elements of the Hamiltonian to the integration over the smaller 
space of its iV eigenvalues. The parameter /3, with values 1,2 or 4 describe 
the symmetry of the ensemble (named orthogonal, unitary and symplectic, 
respectively) . 

If the elements of the random matrix are believed to be statistically in- 
dependent from each other, one obtains the quadratic confinement potential, 
leading to the Gaussian ensembles of random matrices ^3] • The requirement 
of statistical independence is not motivated from first principles. Thus, the 
influence of a particular form for the confinement potential in relationship 
with the Gaussian predictions is an interesting problem [TT] . 

As it is well-known, the joint probability distribution can be mapped onto 
the Gibbs distribution of a classical one- dimensional gas of fictitious particles 
with a pair-wise logarithmic repulsion and a one-particle potential. 

P({ Xi }) =Z- 1 exp[-/3H({x l })] , (3) 
H{{x i }) = - y £hi\x i -x j \+ y £V(x i ) . (4) 

i<j i 

This is usually called the Coulomb gas picture. Let us now present the mean- 
field approximation that was introduced by Dyson [H]. By substituting the 
density of states in the previous Hamiltonian, we obtain a continuous descrip- 
tion of the energy functional H[p] in terms of p(x). The extremum of this 
functional corresponds to an equilibrium of the effective plasma expressed by 
the equation: 

J dx' p{x')ln \x — x'\ = V(x) + c , (5) 

where p is the mean density. This expression is usually known as Wigner 
equation. In jH] , the next correction term was found using a hydrodynamical 
approximation: 

/oo 1/5 
p(x) In | a; — x | dx H — (1 )lnp(V). (6) 
-oo 2 2 
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Higher order terms involve a complicate combination of the density of 
levels with correlation functions [2j. In addition, in jH] it was shown that in 
this approximation the first correction term -the one that appears in (JSJ- is 
of order A^-, while the next term is of order iV -2 . This applies to systems 
with strong confining potential. 

In this paper, we study features of stable random matrix ensembles, in 
analogy with the stable random variables. We provide here a derivation, 
based on simple arguments, of the Cauchy case for the three ensembles {(3 = 
1,2 and 4). In this way, we obtain the so-called Lorentzian ensemble, that 
has been already studied with various points of view [3 IH EH EEj . Further 
work related with stable distributions can be found in |T7j, and in ;3J from 
the point of view of free probability theory 

In section 3, this ensemble allows to provide explicit examples to the 
failure of Dyson equation, in the case of non strongly confining potentials. 
Additionally, we solve several inverse problems using Dyson equation and 
present a systematic comparison between all cases. In particular, we find 
the confining potential that gives a Gaussian density of states in the limit 
N ^ oo. 



2 The stable random matrix ensembles. Deriva- 
tion 

We begin by showing a simple way to obtain the expression for the stable 
ensembles. We will focus on the Cauchy distribution, but all the stable 
distributions can be studied, in principle, with the same method. We look 
for ensembles that satisfy the following: 

N 

M = J2 M i ■ (7) 

i=l 

M and Mj denote the random matrices on the same ensemble. Put it simply, 
the matrix generalization of the stable random variables (see Appendix A). 
Notice the following: 

TV N 

TrM = T T J2M i = Y / TrM; . (8) 

i=l i=l 

By taking the trace of the matrix, we reduce the problem to a problem of 
random variables. In order that this last equality holds, it is evident that 
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the trace itself should be a stable random variable. More clearly, we make a 
mapping between the stable random variables and the stable matrices (the 
matrix counterpart). It is manifest now in which sense we say that we search 
for the matrix generalization of the stable random variables. Note that the 
expression for the Gaussian ensemble satisfy the relation (jHJ), after taking the 
trace on its P(x\, ...,xn). Let us go now into the Cauchy case. We begin by 
writing again the general expression in terms of the eigenvalues: 

N 

P(xx, ...,x N ) = C N \\ \xi - xjf exjp[-^2V(xi)) . (9) 

i<j i=l 

For consistency with the case N = 1, we clearly need the following weight 
function (or confining potential): 

u{x) = exp(-V(x)) = - J— . (10) 

A + x l 

It is also clear that we need some function of iV in our weight function. 
Otherwise, the joint probability distribution would not even be normalizable. 
What we need to impose is: 



/oo /*oo 
••• / dxi---dx^P(xi, xjsf)5(y- 
-oo J-oo 



XN 

- Xi — X2 — ••• — X N) 



with: 



P(x 1 ,...,x N ) = C N ^\x i -x j fl[[j^—^J . (12) 



7T (AA^) 2 + y 2 ' 
(11) 

N_ f i N a(N) 
+ X 2 



From the case iV = 1, we already know that a(l) = 1. We can see that 
the function a(N) should be a linear polynomial in N. Notice that we are 
studying the modification we have to impose to the Cauchy distribution in 
order that the effect introduced by the correlation term does not modify the 
stability property. The Vandermonde determinant introduces the product of 
N(N — l)/2 polynomial factors, while the weight function term introduces 
iV products of polynomial factors. Thus, we need each one of these factors 
to contain a term of order N. So, we have to consider the following function: 

a(N) = a N + ai . (13) 

Then, it is enough to look at the N = 2 case. Needless to say, it can be 
checked that, for example, the case N = 3 implies a zero coefficient for an 
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hypothetical quadratic term, and so on. Actually, in the next section we will 
present an alternative argument that also implies that the function should 
be of the form indicated above. 

Doing the integral for the three ensembles, we finally arrive at: 

2 §(N 2 -N) AT , 1 (l+f(iV-l)) 



p( X1 ,...,x N ) = ** NhN njsi-s/n (ja^j " (14) 

We have arrived at a known and important random matrix ensemble 
[3J EI] • in the physics literature it is named Lorentzian ensemble jl] . The 
normalization constant can be computed with the Selberg integral [T3j . 

2.1. The relationship with the circular ensemble and 
universality. 

It has been shown that the Lorentzian ensemble (|14j) has a density of states 
with same form for any dimension of the matrix and for any (3 4J. Recall 
that the same result holds for the circular ensembles: they have constant 
density on the circle for all N |13j . We know that there are several dif- 
ferent parametrizations for a unitary matrix. A rather convenient specific 
parametrization is S = jt^j (Cayley transform), since it is one on one. 
Notice the effect of any of this transformation on the Cauchy probability 
distribution: 



f(x) = - and x = cot s(^/ 2 )' ( 15 ) 



then: 



■ dx , 1 —1 cos 2 9 1 

" ' = ' d6 ' 7T(1 +COtg(^/2) 2 ) = ' 2C0S 2 ^ ' 71 = 2^' ( ' 

That is, the Cauchy distribution transforms into a uniform probability 
in the circle. In jH , it is shown that under the Cayley Transform, the 
Circular ensemble maps into the Lorentzian ensemble. Consequently, taking 
into account the universality property for the correlation functions of the 
circular ensembles in the limit iV — » oo the Lorentzian ensemble also 
exhibits universal behavior in this limit. That is to say, all the correlation 
functions of the Lorentzian ensemble behave as the ones from the Gaussian 
ensemble in the limit of large matrices [3] . 
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Additionally, notice that if X is a random variable with probability dis- 
tribution function fix) = ,, 2 A 7 . , then the random variable Y = X^ 1 has 

the probability distribution = ^rj^ • Namely, the inverse of a Cauchy 
distribution is a Cauchy distribution. Again, we see that the same property 
holds for the matrix generalization, and in ^H] it is shown that if M is a ma- 
trix in the Lorentzian ensemble, then M _1 is also a matrix in the Lorentzian 
ensemble. This result readily follows by the corresponding change of variables 
in the multidimensional probability distribution function for the eigenvalues. 

Therefore, a matrix in the Lorentzian ensemble shares many properties 
with a diagonal matrix, whose elements are identically distributed Cauchy 
distributions. First of all, the diagonal matrix, after trivial integration over 
N — 1 eigenvalues, also posses a Cauchy density of states (independently of 
N). Secondly, at the level of the density of states it also has the correspon- 
dence with a uniform density. Furthermore, the diagonal matrix is equal to 
its inverse as well, as it readily follows from the property of the Cauchy dis- 
tribution. Needless to say, the correlation properties of the eigenvalues are 
very different. These properties are not shared by a random matrix in gen- 
eral and make the Lorentzian ensemble very useful for physical applications 

El- 

It is also interesting to discuss a little bit more about universality. As it 
is well-known, any strongly confining potential (that is, growing at least as 
V (x) ~ |x| for x — > oo), gives the same Gaussian universality class jE] (see 
also [1J). Here, we deal with a special confining potential, that contains the 
dimension of the matrix N in a explicit way: 

Ve(x) = ( (3N+ 2 2 - p )M\* + x>). (17) 

Note that the dependence with the variable x is given by a very soft confining 
potential. Nevertheless, it is multiplied by a function that is linear in N, 
the dimension of the matrix. This is important, since any random matrix 
universality regime is understood in the limit iV — > oo. Thus, at the end, it is 
the term in N the one that achieves the universal behavior. This clarifies the 
mechanism that makes this ensemble to belong to the Gaussian universality 
class. The consequence is that the ensemble exhibits very different global 
properties (the density of states), while it has the same behavior at the local 
scale (all the correlation functions). 

It is patent that all the arguments above hold for any stable distribution 
(see Appendix A). As an extension, note that this matrix generalization of 



6 



the Levy distributions can be expressed, in general, as: 
P(x u ...,x N ) = C N Y[\xi - xjf n {f a ,i3'(x)) ° \ a + ai = l. (18) 

i<j k=l 

Needless to say, it is consistent that the only stable distribution with finite 
variance, the Gaussian distribution, is the only one that has different form 
for a(N) (more precisely ao = and a± = f3) in the matrix generalization. 

We see that the coefficient a is actually equal to a = @- = -£^. Needless 
to say, as explained above, the value of the coefficient depends both on the 
power of the polynomial in the numerator (the Vandermonde determinant) 
and the one in the denominator (weight function). Thus, this particular 
value seems to be expected for all the Levy distributions. Since we know the 
asymptotic behavior of the distributions for x large (see Appendix A): 

f«A*)~r^ «<2, (is) 

\x\ 

we have the following confining potential for x large: 

V(x) = (^p[ N + ^ In (l + x 1+a ) w (pN + In x , (20) 

and, for N large and x large: 

V(x)=(3N\nx. (21) 

This coincides with the result found in [7. for the Hermitian case ((3 = 2), the 
one considered there. Notice that, at this level of approximation, the con- 
fining potential is the same for any Levy distribution. Since this is precisely 
the universality limit, it is expected that they all share the same universal 
properties in the sense explained above. 

3 The Dyson equation 

We use now the previous results to illustrate the failure of Dyson's mean- 
field approximation in certain cases. This has been recently discussed in the 
literature, and has proved to be an increasingly interesting topic, as soon as 



7 



new type of random matrix ensembles, different from the Gaussian ones, are 
studied. As mentioned, Dyson equation is jH] (see also [2]): 

f°° 18 
V( x ')= p(x)ln\x-x'\dx + -(l-^)lnp(V) . (22) 

We would like to study several nontrivial examples that allow to illustrate 
the situation. In the progress of doing this, we will also solve some interesting 
inverse problem in random matrix theory. 

The Cauchy case allow us for a simple and enlightening example. First, 
we begin by asking which is the confining potential that give us a Cauchy 
density of states in the limit N very large 1 . This reduces to compute: 

T w \ A r°° \n\x — x'\ , A , ,,n 9N N, 9N . , 

V ^ = ~ / '2 dx = T- ln ( A + x = T ln(A + x • (23) 

To begin with, we find that for the Hermitian case (/3 = 2), this result 
gives just half of the right value, since the first correction term is identically 
zero for 3 — 2. Thus, we already see that half of the value would be contained 
in the rest of the following correction terms and we should remember that 
the first one of these corrections is of order N~ 2 in the strongly confining 
paradigm. In fact, with the confining potential obtained in (|23|) , we do not 
even have a normalizable probability distribution function, so it can not even 
define an ensemble. We see, in this way, a clear example for the failure of 
the mean-field approximation. It is interesting to summarize what we obtain 
for the three symmetries: 





Exact 


Mean — field 


3 = 1 


r 2 +1 )ln(l + ^) 


(V)ln(l + x 2 ) 


8 = 2 


Nln(l + x 2 ) 


f ln(l+x 2 ) 


8 = 4 


(2N- l)ln(l + x 2 ) 


(f + |)ln(l + x 2 ) 



(24) 



Notice that the mean-field analysis is still able to give the right analytical 
expression for the part in x of the potential, and again hints that we should 
expect a linear polynomial for the part in N. Nevertheless, it completely 
fails to give the right values for the coefficients of the polynomial. 

The typical use of the Wigner integral is to show, for example, that a 
parabolic confining potential give rise to the Wigner semi-circle law ^3] • This 

1 Note that we already know the answer, is the confining potential in the previous 
section. 
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is standard material, but not completely trivial. Now, the reader might won- 
der what happens with some intermediate cases, like for example a density 
of eigenvalues without compact support but coming from a strongly confin- 
ing potential. A natural choice is a Gaussian distribution for the density of 
states. This is also an interesting problem for other reasons: on one hand, 
since we know in detail the Cauchy case -that turns out to be rather simple- 
then, the Gaussian case naturally comes to mind. On the other hand, a 
Gaussian density of states has a long story in the theory of random matrices 
and it is indeed physical . We have to compute: 

roc „ 

V(x') = / e~ x \n\x-x'\dx . (25) 

This is a non trivial computation, relevant in other fields as well (including 
disordered systems [5] and determinants of random Schrodinger operators 
[12]). For example, it appears in [T2] , where it is left unsolved. We carry out 
the computation in Appendix B. The result is: 

V(x)=x 2 2 F 2 (l,l;^,2;-x 2 ) . (26) 

Note that we need a rather complex expression for the confining potential 
in order to have an ensemble with Gaussian density of states in the limit 
N ^ oo. 

In spite of this complex solution -in great contrast to the simplicity of the 
Cauchy case-, it would be of interest to study it in connection with the two 
body random matrix ensembles. These are old ensembles, that have been 
the subject of considerable revision and renewed interest (see jB] for a recent 
review). They are designed to include correlations among the different matrix 
elements and one of their main features is that they show, in a certain regime, 
a Gaussian density of states. In contrast, we have obtained the confining 
potential that give us this density of states for ensembles with the typical 
symmetries (/3 = 1, 2, 4) of random matrix theory. 

To conclude this section, notice the different results given by the first 
corrective term in Dyson equation depending on the nature of the density of 
states: 

p{x) = Cexp(-a^) -> lnp(x) = InC - ax^ , (27) 
p( x ) = - 1 -> \np(x) = In - - ln(A 2 + x 2 ) . (28) 

"77" A 4- T*^ 77" 
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The constant terms are dropped from the potential since they are reab- 
sorbed in the normalization constant. Observe that the parameter appears 
explicitly multiplying the potential in the stretched exponential case (in this 
case, the parameter is essentially the variance and it is directly related to 
N). In contrast, in the Cauchy case, the parameter 2 drops out, and does 
not enter into the expression of the potential. We see then a big difference 
between the two cases and a complete and systematic study of all the possible 
cases seems worth studying. 

4 Conclusions and Outlook 

We have shown with simple arguments how to arrive to random matrix en- 
sembles that are the matrix generalization, in the sense of random matrix 
theory, of the scalar Levy distributions. In the particular case where an exact 
solution is known (the Cauchy case), an explicit form for the ensemble can 
be obtained. We have arrived to known and well-established results in the 
mathematics and the physics literature j3J QUI El- More important, our 
intuitive method of derivation clearly shows some of the particular features 
of these type of ensembles. It reveals that we need an explicit and, at most, 
linear dependence of the confining potential with the dimension of the ma- 
trix. Namely, the confining potential is so weak that we are forced to this 
dependence. 

Thus, a purely logarithmic term (as exactly appears in the Cauchy case, 
and asymptotically in the others), yields an algebraic expression for the 
weight function, that has to battle in equal conditions with the Vandermonde 
determinant. It is precisely this point what makes this ensemble so special. 
It can also be seen within the Coulomb gas picture, where the Vandermonde 
determinant is a logarithmic two-body repulsive term and the weight func- 
tion is essentially a logarithmic (one-body) confining potential, that needs 
the explicit dependence with N to be able to compensate the repulsive term. 
Then, it is rather natural that we get the same expression for the density of 
eigenvalues for all N, since we have a modified confining potential for each 
N. This is exact for the Lorentzian ensemble and also for the tails of the 
distribution of the eigenvalues for the other Levy cases. 

2 That in the Cauchy case is not the variance, since it is infinite. This is usually referred 
as lack of characteristic scale. 
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The known results on the correlation functions of the Lorentzian ensemble 
are also noteworthy. In addition with the previous arguments, they lead to 
the expectattion that all the Levy ensembles possess the universal behavior. 

We have also seen, in the discussion on mean-field approaches, the differ- 
ences with the Gaussian ensembles and with the strongly confining paradigm 
in general. In particular, the computation of the confining potential that 
gives a Gaussian shape for the density of states in the limit N — > oo is also 
useful to illustrate this. In the Levy cases, due to its power-like behavior we 
need an explicit modification of the potential for each dimension of the ma- 
trix. This leads to a relative simple solution for the density of states that is 
essentially insensitive to the level repulsion (by construction). This is exact 
in the Lorentzian case, and in great contrast with the Gaussian ensembles, 
where its shape goes from a Gaussian for N = 1 to a semicircle when increas- 
ing N. We have also studied the inverse case: after all the influence of level 
repulsion among all the eigenvalues (in the limit N very large), what do we 
need to end up with a Gaussian shape for the density of states ? The solution 
is a very complex expression for the confining potential. This complexity is 
in contrast with the simplicity of the Cauchy case. One interesting remark 
is that a power-like weight function is a natural antagonist for a power-like 
repulsion term. 
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A. The stable random variables 



The fundamental importance of the normal distribution is due to the Cen- 
tral Limit Theorem which is a consequence of the Bernoulli and de Moivre- 
Laplace theorems, and the law of large numbers. 

For our discussion it is important to have clear the following basic results: 
According to Levy, a distribution F is stable if and only if, for the two 
positive constants c\ and c 2 , there exists a positive constant c such that X 
given by: 

CiXx + c 2 X 2 = cX , (29) 

is a random variable following the same distribution F, as the indepen- 
dent, identically distributed random variables X\ and X 2 . Alternatively, if: 

<p(z) = (e iXz ) = / e iXz dF{X) , (30) 

denotes the characteristic function of the distribution F, then F is stable 
if and only if: 

Lp{c x z)ip[c 2 z) = ip(cz) . (31) 

The most general definition can be found in Let X,Xi,X 2 , ■■■,X n be 
iid random variables with a common distribution F. Then F is called stable 
iff and only iff there exist constants c n > and 7 n such that: 

Y n = Y,X t = c n X + ln . (32) 

i 

Then, the characteristic function according to the definition above satisfy the 
functional relation: 

<p n (z) = <p(c n z)e^ z , (33) 
which can be solved exactly, and the result is: 

Proposition 1 

ip(z) = logip(z) = i^z — c \z\ a jl + ifl'yj w ( z i a )| j (34) 

where a, (3', 7, c are constants (7 is any real number, < a < 2, — 1 < (3' < 1, 
and c > 0, and, 
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u(z, a) 



tan if ifa^l 
f log|z| ifa= 1 . 



(35) 



a is called the Levy index or characteristic exponent. The limiting case a = 2 
corresponds to the Gaussian. For — the distribution is symmetric. 7 
translates the distribution, and c is a scaling factor for X. So, these last two 
parameters are not essential and one can disregard them. 



Proposition 2. The asymptotic behavior of a Levy stable distribution 
follows the inverse power-law: 



A 



a,/3' 



\X\ 



|l+« ' 



a < 2 



(36) 



Proposition 3. The analytic form of a stable law is given through the 
Fox function: 



fcuBf (x) = eHH 



x 



(l-e,e) (1-7,7) 

(0,1) (1-7,7) 



for a < 1. 



2:2 



X 



-1,1 

(-€,€ 



-7,7) 
-7,7) 



for a > 1. With the abbreviations e = 1/a and 7 = (a — (3') /2a. 



(37) 



(38) 



Examples: For a = 2, then (3 1 = 0, and the stable density is identical 
to the Gaussian distribution. 

For a = 1 and = 0, the stable density is identical to the Cauchy or 
Lorentz distribution: 



fi,o( x ) 



A 



n(X 2 + x 2 ) 



(39) 



Further examples can be found in El • 



15 



B. Integral computation. 

We need to compute the following integral: 

V(t) = / ln\x-t\e- x dx . (40) 



These type of integrals represent a characteristic polynomial computation 
(the characteristic polynomial associated to a certain continuous density of 
states). Its relevance in different fields of mathematics and physics can be 
appreciated by consulting H2] • 

First, we consider the Taylor series expansion: 

V® = ±^. (41) 
and we make the following change of variables x = y + t, then we have: 

V(t) = l°° ln\y\e- (v+t)2 dy . (42) 

J — oo 



and we consider: 



Art 

( e~ z2 = (-l) k e- z2 H k (z) . (43) 



dz k 

Then the Taylor expansion looks: 

VW = ±t$m>, ,44) 

where V k (0) = X^ln \x\ e~* 2 H k (x)dx. Because of the symmetry of the Her- 
mite polynomials, we arrive at: 



v[t) = j2 1 \. 2k[) t k V 2k (0) = 2 / H 2k (x)e-* 2 Inxdx. (45) 

1,—r, Kl JO 



k=0 

The point is that last integral can be solved and gives: 



We obtain then: 



^(0) = -^-(-l) fe 2 ife r(A;) . (46) 
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Using the recurrence and duplication formulas for the gamma function we 
can arrive at our final result: 

V(t) = (7 + 2/n2 - 2t 2 2 F 2 (1, 1; |, 2; -t 2 )) . (48) 
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